Phonon switching and combined Fano-Rice effect in optical spectra of bilayer graphene 
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Recent infrared measurements of phonon peaks in gated bilayer graphene reveal two striking 
signatures of electron-phonon interaction: an asymmetric Fano lineshape and a giant variation 
of the peak intensity as a function of the applied gate voltage. In this Letter we provide a unified 
theoretical framework which accounts for both these effects and unveils the occurrence of a switching 
mechanism between the symmetric (Eg) and anti-symmetric (Eu) phonon mode as dominant channel 
in the optical response. A complete phase diagram of the optical phonon response is also presented, 
as a function of both the charge density and the bandgap. 
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Single and multi-layer graphenes are among the most 
promising systems for the development of carbon-based 
devices in electronics. Bilayer graphene is of particular 
interest because a controlled tunable gap in the electronic 
spectrum can be there induced by applying one (or more) 
external gate voltages , as observed by transport 3] 
and optical measurements The potential interest 

in application has triggered also an intense research on 
the vibrational properties, that can be used for instance 
for a careful characterization of the number of layers, 
the charge doping and the amount of disorder 
Large part of the work in this context has focused on 
the properties of the in-plane Eg mode at a; w 0.2 eV 
which is present in single as well as in bilayer systems, 
and which can be probed by Raman spectroscopy 

Recently, phonon peaks in the energy range w ~ 0.2 
eV have been reported also in infrared (IR) optical mea- 
surements of bilayer graphene 15, 1^, showing a rather 
different phenomenology with respect to Raman spec- 
troscopy. In particular, the observed IR phonon peak 
presents a strong dependence of the intensity and of the 
Fano-like as ymm etry as a function of the applied gate 
voltage 15|, ll6[. These features have been attributed 



respectively to a charged-phonon effect for the Eu anti- 
symmetric (A) mode 15], or, at n = 0, to the emergence 
of a Fano profile for the Eg symmetric (S) mode, whose 
optical activity can be triggered by the electrostatic po- 
tential difference A between the two carbon planes [l6|. 
However, the possible connection between these two al- 
ternative views is still unclear. 

In this Letter we provide a unified microscopic frame- 
work that allows us to elucidate the relative role of the 
Eu and Eg phonon modes in bilayer graphene, with re- 
gards to the infrared activity and the Fano asymmetry 
of the observed phonon peaks. We present a complete 
phase diagram for the strength of the phonon modes and 
their Fano properties as functions of the chemical poten- 
tial and A, showing that a switching mechanism between 
the dominance of the Eu or Eg mode can be controlled by 
the external gate voltage. Our work permits thus recon- 
ciling within a unique approach the phonon-peak features 



observed by different experimental groups 15, 16j . 

To compute the conductivity of bilayer graphene 
we work in the 4x4 basis of the atomic orbitals 
*k = ('^ik'^ik'^L'^L)' where al^^ and bl^^ operators 
create an electron in the layer i and on the sublattice 
A or B, respectively. In this basis, the Hamiltonian for 
bilayer graphene near the K point reads [l^: Hi^ = 

[hvk ■ !{&) + (A/2)a,(/) + (7/2) [a^ia,) + (Ty(ay)]}, 

where ai and / are 2x2 Pauli matrices and the unit 
matrix respectively, and A{13) = A (E) B- Here v is 
the Fermi velocity for single-layer graphene and 7 is 
the interlayer hopping. The electrostatic potential 
difference A induces a gap in the diagonalized bands 
Ck.n Q, labeled according to Fig. [T]i. We set the 
electric field of the infrared radiation along the x-axis, 
so that the electric current reads j — *kJk^k where 
jk = —evl{ax). The two Eu and Eg in-plane optical 
phonons have degenerate longitudinal and transverse 
polarization at zero momentum q = [l7| . We can 
write the electron-phonon interaction for these modes 



where 



is the dimensionless 
lattice displacement for the i/ —A.S branch at q = 0, 
Vi, = ^k^'^^k is the corresponding electron-phonon 
scattering operator. For a choice of the longitudinal 
polarization along x, one finds Va — igo'zi^'y) and 





Y , ■■>■ Xjj ■■>■ ■>--^XjA 

3 k ^ ^ 
■ =j = Va 



■■>■ 



= Va 



FIG. 1: (a) Scheme of the band structure, (b)-(d) Relevant 
diagrams entering the optical conductivity for A = 0: Dashed, 
solid and wavy lines represent the photon, the electron and 
phonon Green's functions, respectively. Squares and circles 
are the current and the electron-phonon scattering matrices 
j, Va, V2' respectively. 
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Vg — wH^y) fl^i where g is the electron-phonon 
couphng. 

Let us consider first the case of ungapped bilayer 
graphene (A = 0), where only the A mode is opticaUy 
active. Since graphene is a non-polar system, the bare 
dipole induced by the rigid shift of the valence charges 
upon the A lattice distortion is extremely small [15| . 
Nevertheless, an optical phonon response can be still 
mediated by the conduction charges in the presence of 
the electron-phonon interaction, as suggested by M.J. 
Rice in the context of organic and fullerene compounds 
[isl . The total optical conductivity of the system is com- 
puted as cr{uj) = — Xjj(cj)/icj, where Xjji^) = 
is the current-current response function. In the pres- 
ence of electron-phonon interaction, one can identify 
two classes of contributions in Xjji^)^ namely Xjj{^) — 
Xjj(w) + x^j{uj). The former, depicted in Fig. Ht) de- 
scribes electronic excitations (see, for example, Ref. |19j). 
while the latter contains all the diagrams which can be 
spli t in two by cutting one phonon propagator (Fig. [TJ;) 
llSll . We can write thus 



Xjji^) = XjA(w)£'AA(w)xAtj(t^), 



(1) 



where Xja(w) = -(jVa)cj, XAtj(w) = -{Vlj)uj are the 
mixed current-phonon response functions, andZ)AA(w) = 
— ((/)a0a)lj ~ — OJA + iV is the phonon propa- 
gator with frequency oja and linewidth Fa renormalized 
by the phonon self-energy XAtA (Fig- HJl)- R should 
be emphasized that two different response functions, Xj' 
and XjA, enter in the above decomposition of <j{uj). This 
distinction, that has been neglected in the original formu- 
lation 18[, is however crucial, because it implies that the 
allowed particle-hole excitations of the system will con- 
tribute in a different way to the electronic optical back- 
ground, related to xfj^ to the phonon renormalization, 
controlled by Xa+Ai or to the electron-phonon optical re- 
sponse, controlled by Xjh- 

Eq. ([1]) leads to the onset, in the real part of the optical 
conductivity UopC'^) = ~I™.Xjj (i^)/a;, of a phonon peak 
at WA. Indeed, using the relation XAtj(w) = Xjh{^)-, and 
taking the real and imaginary part of each element in Eq. 
dl]), we get 



2 [x;-a(^a)] 



2 ^ 



waEa 

where z = 2[lu — waI/Ea and where 

Xj-a(^a) 



<7A 



(2) 



(3) 



Eq. ([Jlhas exactly the same structure as the Fano for- 
mula |20|. The derivation of Eqs. ©-([3]) shows not 
only that in bilayer graphene the Fano effect originates 
from a correct implementation of the charged-phonon 
Rice theory, but it provides a compelling procedure to 
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FIG. 2: (Color online) Wa, W'a intensities of the _E„ mode as a 
function of the charge concentration n for A = andn = 0, 30 
meV. Also shown are experimental data from Ref. [15| for the 
bare intensity Woxp and the ct)-integrated one VK^'xp. 



evaluate on microscopic grounds the shape and the in- 
tensity of the phonon peak. For instance from Eqs. p|)- 
we can identify the w-integrated peak area as = 
7r[xfA(^A) - x';1{^a)]/^a = (1 - l/gi)[7rxfA(^A)M]. 
Since the sign of W'j^ defends on it can be conve- 
nient, as done in Ref. [15|, to define a "bare" inten- 
sity as Wa = 7rx^^(wA)/wA, which coincides with W'a 
in the limit I^aI — >■ oo, when Eq. ([2]) reduces to a con- 
ventional Lorentzian peak with weight Wa- Note how- 
ever that in the opposite case |(?a| ~ one recovers from 
Eqs. (HI)-® a completely negative Lorentzian peak, of 
intensity — 7r|x"A(wA)P/wA, so that the definition of 'in- 
tensity' of the phonon peak can be ambiguous in the 
presence of the Fano effect. A more convenient quan- 
tity to parameterize the strength of the phonon peak 
is PA = 7r|xjA(wA)P/wA = VFa(1 + l/gi); which is al- 
ways positive and vanishes when the phonon peak is com- 
pletely absent. 

Eqs. (HI)-© can be computed analytically by using the 
non-interacting electron Green's functions. The quantity 
XjA has the typical structure of a particle-hole Lindhard 
response function, with proper coherence factors C"^ 
weighting the contributions of the various excitations be- 
tween the n and m bands. In particular, using the explicit 
matrix expressions of the j and Va operators, one gets 

xM = x]i(^) + x]i(c^)-x,'i(c^)-x|i(c^), (4) 

where xjK^) = 7r^r(w) - ^,T(^)> and 

Here f{x) = l/[exp(x/T) -I- 1] is the Fermi function, 
CfA^k = igevNsNy)^/4:^y (hvk)'^ + 7^ for n,m as in Eq. 
(j4|) and zero otherwise, ^ is the chemical potential, 
Ns = Nv = 2 are the spin and valley degeneracies, and 
the 77 factor takes into account broadening effects due to 
impurities and inhomogeneities. In the clean limit, rj = 0, 
we obtain the analytical expressions valid for T = and 
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ImI < 7: 



Xj-a(^) 



Alln 



{l + u)il-u + 2w) 



{1- u){l + u + 2'w) 



Auw 

(1 - U2) 



(6) 



x'-aH = 7rA[e{\u\- 1)9(1 -\u\ + 2w)+2wS{\u\ 
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where A = gej/Anhv, u — fiw/^ and w = |/i|/7. Similar 
analytical expressions can be obtained for /i > 7. 

To compute the optical conductivity ([2]) we evaluate 
the complex function XjA('^) &t the phonon resonance 
wa ~ 0.2 eV, using standard values for fm — 6.74 
eVA and 7 — 0.39 eV. From the value a ~ 6.4 cV 
[2l| of the deformation potential we get g — 0.27 
eV, corresponding to the dimensionless parameter A — 
{■\/^/TT)g'^ /{hv /a)'^ = 6 X 10~^ [l7[, in agreement with 
the experimental estimates given in Refs. [l^, 15 1. The 
resulting spectral weight Wa as evaluated from Eq. ^ 
and Eqs. (U)-® is shown in Fig. [5] for 77 = and 77 = 30 
meV, alon g wi th the experimental data Wcxp, W^oxp taken 
from Ref. [15|. Note that the only possible electron- hole 
excitations at the phonon energy uj^ (< 7), namely the 
2-3 multiband transitions, are not allowed in Eq. Q. 
Thus in the clean limit Xj'aI'^a) ~ 0, so that no Fano 
effect (gA ~ —00) is found and W'^ = W^, while for 
77 ^ a small contribution to ^"^(wa) from the other in- 
terband transitions gives rise to a weak Fano asymmetry, 
and hence ^ Wa , as shown in Fig. [2] Both the mag- 
nitude of Wa , which is proportional to the dimensionless 
coupling A, and its doping dependence are in excellent 
agreement with Woxp, pointing out that the Eu mode is 
the main responsible for the phonon infrared intensity 
reported in Ref. [l5| at large n. On the other hand, even 
in the presence of a large 77, the above calculations do not 
account for the negative integrated weight Wg'^p observed 
around ri « (Fig. [2]), which has been attributed in Ref. 
[lit to the onset of the Eg mode in the presence of a finite 
potential asymmetry A. These observations suggest thus 
that different phonon modes can be optically relevant in 
different regions of the phase space. 

Our theoretical framework allows us to investigate the 
possibility of such phonon switching by taking into ac- 
count explicitly the role of A. When A 7^ two effects 
must be taken into account: (i) the phonon eigenmodes 
of the systems do not correspond any more to Eu and 
Eg, even though the new phonon eigenfrcqencies uj± fol- 
low closely the doping dependence of the wa,s of the un- 
coupled modes [22, |23| . As a consequence the phonon 
propagator Daa can develop at large A a second peak 
with weaker intensity at approximately the frequency ujs 
of the Eg mode. Note however that according to Eq. 
([1]) the infrared activity of the Daa phonon propagator 
is still ruled by the strength pA- Since, as we shall see 
below, PA vanishes for n also when A 7^ 0, the ap- 
pearance of the IR phonon structures at rt ~ 1^, 16 1 
cannot be accounted for by the mixing of the modes; 
(m) in addition to the previous effect, the presence of 



a finite A leads also to a finite mixed response function 
Xjs(w) = "{jVs)uj 7^ 0. Eq. H]) must be thus generalized 
as: 



XjA{uj)DAA{i^)XA^ji^) 

-t-Xjs(w)£'ss(w)xstj(w) 

+ [XjA('^)-DAs(t^)Xstj('^) + h.c] , (8) 



where Xstj(w) = Xjsi^^), and Di,^,{uj) = -{(t)„(l>u')oj are 
the phonon propagators, calculated including the hy- 
bridized phonon self-energy XAS for finite A 22, [2^. Eq. 
([8]) shows how, due to the Xjs response triggered by the 
finite A, a direct coupling channel to the symmetric Eg 
mode vibrations (the Z3ss phonon propagator) appears in 
the optical spectroscopy. Similar expressions as Eqs. ([4])- 
(O can be derived for XjS j where however the coefficients 
C"g"j., as well as C"j^^, have a more complex structure 
for A 7^ 0. 

To elucidate the competition between the different op- 
tical channels in Eq. ([8]) we compute explicitly XjA and 
XjS for generic A and fi, giving a complete phase dia- 
gram that can be explored in double-gated samples. To 
parameterize the strengths of the relative channels we 
use the quantities Pi, — 7r|xji/p/a;av for v — A,S and 
PAS — "'IXjAXjsl/'^av for the mixed channel, where Wav 
is the average frequency of the two poles of the phonon 
propagators [ 22 . t23| . For a direct comparison with exper- 
imental data, we set T ~ 10 K and r/ = 30 meV, which 
is halfway between rj = 18 meV reported in Ref. Q and 
77 = 40 meV in Ref. [l^. The results for the phonon 
strengths and Fano factors are plotted as a color-map in 
the A-/i space in Fig. [3^-f where also the A-fi location of 
the available IR data 15, 16j are shown. Fig. [3^,b shows 
that the main parameter tuning the strength of the i?„ 
mode is the charge doping whereas the Eg mode is ac- 
tive mainly around the neutrality point, with a strength 
that is tuned by the asymmetry gap A. The strength 
of the mixed optical structure [third line in Eq. ([5])], 
shown in Fig. [31;, satisfies pas ^ Pa,Ps in a large part 
of the phase diagram so that only the A and S channels 
[first two lines in Eq. ([U] are active. The relative inten- 
sity Ras = {pa - Ps)/{pa+Ps) is summarized in Fig. 
[2Ji where Ras ~ 1 corresponds a dominant £"„ mode, 
while Ras ~ — 1 gives a dominant Eg resonance. Note 
that the region Ras ~ where the two strengths of the 
optical structures have similar magnitude is quite nar- 
row and difficult to resolve experimentally. In Fig. [3^,f 
we plot also the relative Fano factor q^, for the A and S 
modes. As one can see, qa is essentially doping indepen- 
dent while qs shows a sizeable dependence as function of 
/X. These behaviors can be understood considering that 
the structure of Eq. ([4|) for XjA is still valid for A 7^ 0, 
so that the low-energy 2 — 3 interband transitions are 
missing and the weak Fano asymmetry of the A mode 
is only due to a finite broadening due to 77 on the re- 
maining transitions. On the contrary all the interband 
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FIG. 3: (color online) (a)-(c) strength of the different 
electron-phonon contributions to the optical conductivity. 
The above color scale is in units of 10^ Jl"^ cm~^. (d) Relative 
phonon weight -Ras. (e)-(f) Fano asymmetry factors gA and 
gs. Also shown in these plots is the A-/.1 location of exper- 
imental data: [l^ (•), (■). The corresponding optical 
spectra are shown in panels (g) (/i = and increasing A from 
A to E) and (h) (increasing doping from 1 to 21), respec- 
tively. Curves are displaced vertically for clarity, (i) different 
electron-phonon contributions Wl (open symbols) to total cj- 
integrated spectral weight (solid line) calculated for A-/i of 
Ref [l5[. (j) corresponding Fano factors g^ from Eq. ([3)1. 
The size of the open symbols in panels (i)-(j) is proportional 
to p^. Crosses are experimental data from Ref. [151] . 



transitions contribute to XjS i including in particular the 
low-energy 2 — 3 interband transitions which overlap with 
the phonon frequency (for 2|^| < ws), accounting for the 
sizable dependence of gs as function of /i. 

The comparison of the present results with the A-/i 
location of the experimental available data provides an 
important route to check our theoretical predictions. In 
Fig. Efe-h we show the optical conductivity o'^^{iJ) using 
Eq. dS]) for the [i, A values corresponding to the exper- 
[II 



imental data of Refs. [15|, |16|. The phonon propagators 
Dyiii are computed takin g ac count the self-energy hy- 
bridization due to A ^ [22, l23[ , and convoluted with 
the experimental resolution of 10 cm~^ [l5l|. In Fig. [3^ we 



plot the spectra for /i = 0, showing the evolution of the 
optical intensity and of the Fano asymmetry for increas- 
ing A 



These features are related uniquely to the Eg S 



mode [16|, consistently with panels (a) and (b). On the 
contrary the spectra in Fig. |3h, evaluated at the dop- 
ing levels and electrostatic potentials of Ref. [l^, are 
expected to show a continuous switching between the 
A and S modes, depending on the values of /i and A. 
To investigate deeper this issue, we plot in Fig. O the 
theoretical spectral weights Wl for both the and Eg 
modes evaluated on the set of experimental data A-/x of 
Refs. [l^. The size of the symbols is proportional to the 
peak strength y^. The switch between the E^ phonon 
peak and the E^ one in the regime n € [— 1 : 3] x 10^^ 
cm~^ is evident and is reflected in a total spectral weight 
W!y = W^+Ws which becomes negative in such n-region, 
in good agreement with the experimental w-integrated 
weight Wg'xp. A similar phonon switch is also evident 
from the analysis of the Fano asymmetry, reported in 
Fig. 13] • Also here one can distinguish the crossover be- 
tween a constant q behavior at large n, which we can 
attribute to the Eu mode, and a drop of |g| at small n in 
the region where the Eg mode becomes dominant. The 
experimental fit, done with a single mode Fano formula, 
presents a similar trend. 

We stress that the switching between the Eu and Eg 
modes in the optical conductivity discussed here is not 
related to the possible appearance at large A of a double- 
peak structure in each of the -Daa or Z?ss propagators 
[13, 14, 22, 23]. Even though the two effects could be 
present simultaneously at large gap values, the optical 
data available so far are outside the region where the 
two peaks can be resolved. 

In conclusion, in this paper we presented a complete 
theoretical description of the phonon resonance in bilayer 
graphene, that accounts for both the intensity and the 
Fano-asymmetry variations as functions of the density 
and the gap. We also showed that an optical switch- 
ing from E^ to Eg can be induced in a controlled way, 
providing a full understanding of the experimental data. 

The derived phase diagram for the optical properties 
offers a roadmap for the characterization of graphenic 
systems. For instance the measurement of the IR in- 
tensity of the phonon peak can provide a useful tool to 
determine the doping level in contact-free samples. 
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